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ABSTRACT 

The problem of transporting recruits from Navy Recruiting Stations 
to Recruit Training Centers in such a way as to minimize total cost is 
subjected to solution by standard linear programming methods, A sensi- 
tivity analysis of the variation of total cost for different modes of 
transportation, as influenced by the variation of the proportion of 
recruits going to each of the two Recruit Training Centers, is carried 
out. A probabilistic model of the occurrence of a given number of 
enlistments, up to a fixed quota, at each Recruiting Station is formu- 
lated and a simulation of the transportation problem arising is carried 
out. A rationale is developed for the derivation of a realistic 
estimator for planning, and some of the difficulties of determining 
a best estimator are discussed. A computer program for approximating 
the expected loss for any estimator is presented. Someother types of 


approach to the problem are indicated. 
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he Introduction, 

In today's Navy considerable emphasis is being put on cost re- 
ductions of all forms. One important example of an area in which 
Significant cost reductions might be made is that of the transportation 
of recruits from Navy Recruiting Stations (RS's) to Recruit Training 
Centers (RTC's). The large number of such "shipments" made indicate 
that close study of the costs are needed. Accordingly, the Bureau of 
Naval Personnel, in their letter PERS B-61 KAD of 29 May 1963, re- 
quested investigation of methods to reduce these costs. This paper is 
presented as an illustration of a practical application of the basic 
linear programming method of solving transportation problems as applied 
to naval recruits from the time of their enlistment until their arrival 
at the appropriate RTC. 

There are thirty-nine main Naval Recruiting Stations and nineteen 
class "A" substations, a total of fifty-eight recruit sources, Each RS 
is assigned, as a quota, a certain percentage of total number of recruits 
to be included by the Navy Recruiting System for the year. To receive 
the recruits there are two RTC's each of which is assigned to receive 
half the total number of recruits enlisted. Recruits can, in principle, 
be sent from any RS to either RTC by any of three modes of transporta- 
tion: air, rail, or bus, all of which cost different amounts. In 
addition, there are differences in cost which are a result of obtaining 
group rates. In order to analyze the cost of transportation more than 
superficially, it is necessary to include costs other than those of 
tickets, e.g., meals, miscellaneous (limousine service, etc.) and pay 
of recruits during transit. The desired solution of the problem 


should give the optimum number of recruits to be sent to each RTC from 





each RS and should determine the modes by which they shouldbe sent. 
In addition, since fare rates, cost of living, pay rates, RS quotas, 
and RTC quotas may all vary, it is necessary to determine the effects 
of all of these changes on the optimum transportation allocations, 
Determinations of these effects are known as sensitivity analyses. 

This problem, with slight modifications, may be formulated as the 
linear programming transportation problem. This formulation of the 
problem is that of the shipment of a homogeneous preduct from a stated 
number of sources to a stated (pessibly different) number of receiving 
points (destinations): the amount to be shipped from each source and 
the amount to be shipped to each destination is predetermined; however, 
the amount going from a specific source to a specific destination is 
not. The cost of shipment of one item from any source to any destina- 
tion is a necessary input. The solution determines the amounts to be 
shipped from each source to each destination in such a manner as to 
minimize total cost, and gives the minimum total cost. 

In this paper we have made use of a previously written FORTRAN 
program of the simplex algorithm to solve (on a CDC 1604 computer) the 
basic transportation problem as posed by the Pureau of Naval Personnel. 
In addition we have investigated the changes in the solution introduced 
by varying certain of the constraints (e.g., the restriction that 50% 
of the men go to each Recruit Training Center). Further, an attempt has 
been made to come to grips with the more realistic problem in which the 
number of recruits desired is specified, but the number which will 
actually enlist can be described only probabilistically. 

Since some of the data were unavailable at the time this paper 


was written, those values were supplied by educated guess. Further, 
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the actual number of men in the quotas were not known; the percentages 
of the total were used instead. Thus, this paper is a demonstration 
of a method which becomes a true solution (within the limits of its 


assumptions) only when the true values are supplied. 





an Problem Formulation. 

In this (classical) statement of the linear programming transpor- 
tation problem the following information is necessary as input: 

It is desired to ship a homogeneous product from each of a given 
number of destinations in such a manner that each destination receives 
an exact desired amount (which can be different for each destination). 
Further, this shipment is to be carried out in such a manner that each 
origin sends an exact desired amount which can be different for each 
origin. The total amount that is desired to send is equal to the total 
amount that is desired to receive. No shipment takes place from destina- 
tions to origins. The problem is to be solved in such a way as to mini- 
mize the total cost of all this transportation. In other problems some 
of the above equalities may be inequalities, but in our problem equal- 
ities are adequate. 

To solve our problem realistically, account must be taken of the 
fact that seldom will money savings be the only governing criterion in 
any problem involving people. Hence, it will pe necessary to provide 
solutions to the transportation problem under consideration, which for 
example, might be to prevent trip length from being excessive. Suppose 
a decision were made (for reasons of morale, say ) that no recruit should 
spend more than thirty-six hours in travel. Obviously for some RS's 
the cheapest mode of transportation, which is also the slowest, would 
no longer be usable, and some costs would be changed. This would raise 
the total cost of transporting recruits and might well change the pattern 
of shipments, i.e., how many men are sent between which sources and 
destinations, for minimizing total cost using the revised costs. Another 


variation introduced into the problem is that caused by the availability 





in some areas of reduced rates when drafts of more than a certain 
number of men are sent at one time. A solution of the problem using 
the assumption that advantage is taken of all such available group 
rates can be justified by assuming that enlistees can be told to report 
for transportation a long enough time after first enlisting to ensure 
that the station has enough recruits to obtain group rates, It should 
be noted that various standard queueing theory models of this situation 
might be used to select appropriate delays. 

The input used in the problems solved for this paper were obtained 
as follows: 

Costs of air, train, and bus tickets from each RS to each RTC 
(San Diego and Great Lakes) were requested from each RS. In addition, 
the RS's were asked to supply meal costs, miscellaneous costs, and time 
enroute for each mode. Not all requests were answered, so missing values 
were supplied by estimates, In the data supplied, the numbers of men 
required at the destinations and available at the origins were stated 
in percents, the information having been supplied by the Bureau of Naval 
Personnel. The letter prompting this study specified that each Recruit 
Training Center should receive the same number of recruits, giving the 
effective determination of that input as fifty percent. 

Once having solved the basic problem for the various conditions of 
restriction on cost (maximum cost, entnee single fare, minimum group 
fare, minimum with time limit, etc.), it would be of interest to inves- 
tigate the effect on the optimal solution when changes are made in the 
system, Since in the transportation problem, the addition of new 
variables (RS's or RTC's) requires not a variation of the same problem, 


but the solution of an entire new problem, such additions will not be 





investigated. The changes in the system which appear to lead to fruit- 
ful investigations are changes in the costs, in the proportion of re- 
cruits assigned to each RTC, and in the numbers available at each RS, 
For a simplified mathematical statement of the foregoing formulation, 


see Appendix A. 





3. Sensitivity Analysis Using Cost Inputs as Parameters. 
3.1 The following assumptions are used for this analysis: 

First = The statistics presented in the tables are NOT OFFICIAL 
NAVY statistics and the conclusions from this analysis are strictly the 
personal opinions of the authors. 

Second ~ It was stated in the introduction of this paper that there 
are thirty-nine main RS's; however, this analysis was conducted using 
forty since we believed at the time that Baltimore was a main Recruiting 
Station. This analysis was completed prior to the discovery of this 
oversight. This oversight has no effect on the conclusions or the 
validity of the described method for solving the problem, Hence, for 
this analysis, recruits are sent from forty RS to two RTC's in accordance 
with the travel rates or cost coefficients listed in Table I, page 25. 
These rates are constantly changing; therefore, this table must be “up- 
dated" to be applicable to a particular time, Computed in these rates 
are subsistence which includes meals, berthing, limousine service, etc., 
and pay during travel which amounts to twelve cents per hour. 

Third - The overall annual Navy recruiting quota is distributed 
among the forty Recruiting Stations as listed in Table II, page 26. 
These percentages were computed by multiplying the annual district quota 
by the individual Recruiting Station's annual quota. For example: 

The first recruiting district is responsible for 16.52% of the overall 
annual Navy quota, Of this percent, Boston is responsible for 27.00% 
within the first district; therefore, Boston is responsible for 27.00 
times 16.52 or 4.46% of the overall annual Navy quota. These percentage 
assignments are redistributed slightly each year; therefore, this table 


also should be "up-dated" to be applicable to a particular year. 
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Furthermore, some of the RS's have class “A" substations which further 
breaks the percentage downs however, for the problem discussed in this 
paper only the forty main recruiting stations were considered as origins 
because the percentage breakdowns to the class "A" substations was not 
available. Again, this does not change the basic method of solving the 
problem, 

Fourth - It is emphasized most emphaticaily that this analysis is 
based only on costs incurred during transportation and ignores any con-= 
Sideration of costs incurred before arrival of candidates at RS's (isec. 
Recruiting Stations operating costs) or after arrival at RTC's \i.e., 
costs of operating the RTC's). 

The sensitivity analysis was conducted by using a fixed set of 
cost input parameters and varying the percent recruit input to RTC 
Great Lakes from 0% to 100%. Since the total of 100% recruits is sent 
to RTC Great Lakes and RTC San Diego, it is evident that RTC San Diego 
must receive the complement of those recruits sent to RTC Great Lakes. 
When this approach was used a particular sequence of recruiting stations 
and corresponding percentages resulted; these percentages were the levels 
of input to RTC Great Lakes at which the corresponding RS's commenced 
sending recruits to RTC Great Lakes. On each of the following four pages 
the parameter of cost input is defined. Below each is tabulated the 


applicable sequence of RS's,. 





Max Cost Inputs: 


Max Cost Inputs are defined as the maximum 


transportation costs, regardless of mode, between a RS and the two 


Recruit Training Centers. 


on the costs, 


These inputs were used to put an upper bound 


The use of these cost inputs resulted in transporting 


the recruits via all three modes; however, each recruit traveled via 


only a single mode. 


Sequence of RS's and the Corresponding Percentage at 


Which Each RS Commences to Send Reeruits to RTC Great Lakes 


New York City 
Milwaukee 
Chicago 

Cleve land 
Indianapolis 
Baltimore 
Albany 
Pittsburgh 
Philadelphia 
Minneapolis 
Ashland, Ky 
Buffalo 
Washington, DC 


Columbia, S.C. 


1% 

8h 
10% 
14% 
17% 
19% 
20% 
23% 
26% 
31% 
34% 
36% 
38% 
39% 


Detroit 
Louisville, Ky 


Columbus 


Jacksonville ,Fla 


Richmond, Va 
Boston 
Nashville 
Cincinnati 
st. Louis 
Omaha 

Little Rock 
Macon, Ga 


Birmingham 


10% 
LA i 
L6% 
L7h 
50% 
51% 
5 5b 
58% 
5% 
62% 
64% 
65% 
66% 


Raleigh 

Des Moines 
New Orleans 
Kansas City 
Oklahoma City 
Houston 
Dallas 

Denver 
Albuquerque 
Seattle 
Portland 

San Francisco 


Los Angeles 


68% 
70% 
72% 
73% 
76% 
77% 
79% 
81% 
83% 
85% 
87% 
90% 
9% - 100% 





Air Cost 


puts: 


transport all recruits only by air. 


Air Cost Inputs are defined as the cost to 


These inputs were used primarily 


to arrive at a cost figure when minimum time enroute is the primary 


consideration. 


For Chicago, where no air transportation is available 


to RTC Great Lakes, the train fare was used under the assumption that 


time spend enroute was less than by bus. 


Sequence of RS's and the Corresponding Percentage at 


Which Each RS Commences to Send Recruits to RTC Great Lakes 


Albany 
Louisville 
Cleveland 

New York City 
Boston 
Ashland 
Richmond 
Philadelphia 
Detroit 
Milwaukee 
Washington, DC 
Chicago 
Pittsburgh 


Buffalo 


1% 

3% 

iy 

1% 
15% 
19% 
21% 
22% 
27% 
31% 
33% 
34% 
38% 
42% 


Indianapolis 
Baltimore 
Columbus 
Nashville 
Cincinnati 
Raleigh 
Columbia 
Jacksonville 
St. Louis 
Macon 
Birmingham 
Des Moines 


Minneapolis 


10 


LL to 
45% 
47% 
48% 
50% 
52h 
53% 
55% 
57h 
60% 
61% 
63% 
65% 


New Orleans 
Omaha 

Kansas City 
Little Rock 
Oklahoma City 
Dallas 
Houston 
Denver 
Albuquerque 
Seattle 
Portland 

San Francisco 


Los Angeles 


69% 
70% 
72% 
7L% 
76% 
77% 
79% 
81% 
83% 
85% 
87% 
90% 
94% = 100% 





MinT Cost Inputs: 


costs allowing for a maximum time enroute of thirty-six hours. 


MinT Cost Inputs are defined as the minimum 


Thirty- 


six hours is an arbitrary value assumed to be realistic by the authors 


and any other value might be used to obtain comparable results. 


With 


these cost inputs, as with the Max cost inputs, the transporting of 


recruits via all three modes resulted, 


Sequence of RS's and the Corresponding Percentage at 


Which Each RS Commences to Send Recruits to RTC Great Lakes 


Ashland, Ky 


Boston, Mass 


Louisville, Ky 


Richmond 
Philadelphia 
New York City 
Albany 
Cleveland 
Detroit 
Baltimore 

Mi lwaukee 
Columbus 
Nashville 


Pittsburgh 


1% 
2h 
6% 
8% 
% 
14% 
21% 
24% 
27% 
31% 
33% 
35% 
36% 
38% 


Washington, DC 


Macon 
Indianapolis 
Buffalo 
Chicago 
Burmingham 
Cincinnati 
Jacksonvi lie 
St. Louis 
Raleigh 
coidnene 
Des Moines 


Minneapolis 


11 


41% 
42% 
43% 
45% 
47% 
52% 
54% 
55% 
58% 
61% 
62% 
63% 
65% 


Little Rock 
New Orleans 
Omaha 

Kansas City 
Denver 
Oklahoma City 
Dallas 
Houston 
Albuquerque 
Portland 

San Francisco 
Seattle 


Los Angeles 


69% 
70% 
71% 
73% 
76% 
71h 
79% 
81% 
83% 
85% 
87% 
92h 
QL% - 100% 





Bus Cost Inputs: 


transport all recruits by bus. 


with these inputs. 


Bus Cost Inputs are defined as the costs to 


No maximum time enroute was considered 


We realize that it is not feasible to order a 


recruit to ride a bus for seventy-elght hours from Portland, Maine to 


RTC San Diego, however, these inputs also correspond to the minimum 


costs available and thus, put a lower bound on the costs. 


Sequence of RS's and the Corresponding Percentage at 


Which Each RS Commences to Send Recruits to RTC Great Lakes 


New York City 
Philadelphia 
Mi lwaukee 
Chicago 
Pittsburgh 
Baltimore 
Indianapolis 
Detroit 

Cleve land 
Minneapolis 
Washington, DC 
Richmond 
Cincinnati 


Boston 


1% 

8% 
13% 
15% 
19% 
22% 
2L% 
25% 
30% 
33% 
36% 
37% 
38% 
10% 


Columbus 
Ashland 
Louisville 
Macon 
Buffalo 
Columbia 
St. Louis 
Raleigh 
Albany 

Des Moines 
Nashville 
Jacksonville 


Little Rock 


IZ 


LAL f 
46% 
47h 
4% 
50% 
52% 
5A% 
56% 
58h 
60% 
62% 
64% 
67% 


Birmingham 
Omaha 

New Orleans 
Kansas City 
Oklahoma City 
Dallas 
Houston 
Denver 
Albuquerque 
Portland 
Seattle 

San Francisco 


Los Angeles 


68% 
70% 
72% 
73% 
76% 
77% 
79% 
81% 
83% 
85% 
87% 
90% 
OL% - 100% 





3.2 A graph of percent recruit input to RTC Great Lakes vs. trans- 
portation costs per one-hundred recruits is given in Figure I, page 14. 
This graph illustrates how the variation of recruit input between the 
two RTC's affects transportation costs, The curves on the graph were 
plotted by using the computed data found in Tables III and IV, pages 


27 and 28 respectively. 


3.3 A FORTRAN listing of the linear program used for the previous 
sensitivity analysis is given in Appendix E, This listing is included 
as an aid for those who desire to pursue this transportation problem 
further; a detailed explanation of this program is not considered 
essential to the results of this paper. A format of the computer 
print-out of the above program for a 5% interval, starting at 30% 
recruit input to RTC Great Lakes, and using 1% increments and Air 

cost inputs is given in Tables VI-A through VI-F, pages 30 to 35, 

Table V, page 29, presents a number code which correlates the P-vectors 
in Tables VI-A through VI-F, pages 30 to 35, to their respective RS's 


and RTC's, 
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4. Statistical Variation. 

In the previous formulation of our problem, the basic assumption 
was that the inputs to the problem were known perfectly, but it is 
Quite likely that this will not be the case. Indeed, it is apparent 
that there is a positive probability that the number of men enlisted 
at any RS will not be equal to that RS'’s assigned quota. In other words, 
the number of men to be transported from any RS may be regarded as a 
random variable. Regarding this random variable, /S, as a parameter 
to the problem, we note that Garvin 3 | points out that for any 
parameter, B , whose probability density (or mass) function is known, 


the minimum of the objective function, z,, obeys the relation 


E[z, (B)| Sz, (e[B)) (4.1) 


where $ 
|B | is the vector of expected values of the B's. 
zZ, (B) = min z(b), the minimum objective function expressed as a 
function of B,. 
The lack of equality above makes the problem of saying something 
meaningful about z more difficult. 

Even if we could easily find the value in the left side of equa- 
tion (4.1), we would have an expression only for an ideal situation 
which cannot exist in the real world. In fact what we wish to know is 
what policy we should choose in order to approach this ideal more 
closely than by the use of any other policy. One approach to such 

2 


problems is to find an estimator for the parameter B, and then to use 


this estimator, EP to define our policy. 


AS 





In attempting to find the optimum estimator, B, we have made 
certain assumptions concerning the uncertain (i.e., statistically 
varying) value of B, which is a vecter representing the number of en- 
listments at each RS plus the number of enlistees going to one of the 
RTC's. These assumptions and the models arising therefrom are dis- 
cussed fully in Appendix B. 

It was found (see Appendix B) that there is no simple analytical 
solution to the problem of determining a best estimator, The problem 
of even evaluating any given estimator by analytical means is so great 
as to be impractical to solve; so the authors have taken the approach 
of using the same assumptions in devising a Monte Carlo simulation. 

As discussed more fully in Appendix B, it is assumed that the number 

of people arriving at a recruiting station follow the Poisson probabil- 

ity law and that the number of those who actually enlist is governed 

by the binomial probability law. In our simulation, this random amount 
is determined for each RS and the problem solved to find the least cost 

for that set of inputs. The program is designed to solve any desired 

number of such subproblems in order to find an approximation to the 

value of the left side of equation (4.1). 

The next assumption we make is that there is, conceptually, some 
loss in not solving the sub-problems this way, and that by following 
a fixed policy we will have to incur this loss at least part of the 
time. By the use of an estimator we can find a policy. The policy 
is the transportation pattern developed by solution to a sub-problem 
in which the estimator that we choose is used as the input. When this 
policy is applied to the sample of sub-problems we obtain, for each 


sub=problem, a difference in cost which will be either zero or more 








and which is the loss for that particular sub=problem, In our Monte 
Carlo simulation we have evaluated only one estimator, the average B, 
Of course, given any other estimator, we could evaluate it also, 
Having found the losses we can find the average loss, which is termed 
sample risk, We can also find the sample deviation of loss, thereby 
gaining some idea of how much any particular loss may vary from this 
average, With larger and larger numbers of sub-problems we can more 
closely approach the "true™ values of these statistics. The results 
of the Monte Carlo simulation are discussed in section 5.3. 

Another possible estimator which might be used under certain circum= 


stances is discussed in Appendix D, 








5. Results and Discussion. 

5.1 It is interesting to observe that the sequence of the RS's 
changed when the cost inputs were changed for some reason, such as 
changing the mode of travel (air to bus), or changing the criterion 
for choosing a particular mode of travel (limited time enroute). 

Using Max Cost Inputs, it is noted that when the input to RTC 
Great Lakes is changed to 1%, New York is the first RS to send recruits 
to RTC Great Lakes, and it sends 1% of its quota. New York's total 
quota is 7.35%; therefore, as the input to RTC Great Lakes is increased 
by 1% increments, New York increases its quota by 1% increments until 
its total quota is sent there, One explanation for this could be that 
New York is responsible for the largest cverall Navy quota. However, 
referring to Table I, page 25, it is observed that the DIFFERENCE in 
Max costs for New York is $211.02 = $44.96 = $166.06, This difference 
is larger than for any other city; therefore, it would seem that this 
combination would give the greater savings. 

As the recruit input to RTC Great Lakes continues to change by 1% 
increments it is noted that when the input reaches 8%, Milwaukee is the 
next RS to change its distribution and commence sending recruits to 
RTC Great Lakes. Why should Milwaukee send recruits to RTC Great Lakes 
before Chicago? Once again it is observed that the DIFFERENCE in Max 
costs for Milwaukee is $152.99 = $10.09 = $142.90, compared to $143.39 
- $1.37 = $142.02, for Chicago. It is interesting to note that the 
input to RTC Great Lakes must reach 10% before Chicago starts sending 
recruits to Great Lakes. Although Chicago is the closest recruiting 
station to Great Lakes, this does not impiy that the minimum distance 


traveled is going to give maximum savings, 





Again referring to Table I, page 25, only this time using the 
Air or Bus costs, it is observed that the first RS to commence sending 
recruits to both RTC's is not New York, but that city that has the 
greatest DIFFERENCE in the appropriate costs. When that RS has com- 
menced sending its total quota to RTC Great Lakes, the next city with 
the greatest DIFFERENCE in costs commences to send: recruits to both 
RTC's. This method continues until all RS's are sending their total 
quota to RTC Great Lakes. Therefore, when the problem is linear and 
involves two destinations, as in this problem, we see that it is the 
DIFFERENCE in costs of the applicable mode of travel to the two des- 


tinations that is the controlling factor. 


5.2 The graph in Figure 1, page 14, illustrates that there is a de- 
crease in travel costs up to an 82% input of recruits to RTC Great 
Lakes and then the costs begin to increase. In the case of Max and 
MinT cost inputs, where it is recalled that mixed travel modes were 
used, the minimum occurred at roughly the same percentage. Thus, the 
decrease occurs regardless of the mode of transportation. 
Why should the minimum costs occur at such a high percentage 
input to RTC Great Lakes? It can be conjectured that one possible 
explanation might be: 
1. It costs less to transport recruits from all parts of the 
country to RTC Great Lakes than to RTC San Diego because RTC 
Great Lakes is geographically more centrally located. 

2. The more populous areas are responsible for a larger overall 
-annual recruiting quota and the resonant of these areas are 


located nearer RTC Great Lakes than RTC San Diego. 
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We would like to emphasize at this point that we have discussed 
only that the transportation costs to RTC Great Lakes are less than 
those to RTC San Diego up to a certain percentage of recruit input. 
We have made no cost analyses of sending recruits through the two 


Recruit Training Centers. 


5.3 In testing the Monte Carlo approach to this problem, first a 
sample of fifty randomly generated B-vectors were analyzed. The inputs 
to the random vector generator sub-program were, 

1) RS quotas expressed as integer numbers, 

2) the average number of interviews held per day at each RS, and 

3) the individual probability of success (obtaining an enlistment) 

per interview for each RS. 

Values of the last two parameters were quite Saat, chosen, the | 
average interview number being made approximately a fourth more than 
the annual quota divided by 365. The probability of success at an 
interview was chosen arbitrarily over a range from 55% to 80%. A table 
of these inputs is given in Table VII-B. The vectors resulting had | 
entries of a fair number of stations (about half in most cases) which 
fulfilled their quota. A typical vector resulting is shown in Tatle VII-C. 
The other percent-of-quotas-enlisted ranged down to 64%, but the total 
quota was filled to a level never less than 93% aa never as great as 
95% in the fifty samples. Further, the Pein pattern of transporta= 
tion was the same for each random er for these arbitrary parameter 
values. Typical results of the solution of one vector are given in 
Table VII-D, In short, the solution was quite insensitive to changes 
generated in a small sample by these parameter values, and the sample 
aevietion of the random variable 6 for the RTC was small. Because 
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of this result, it was decided that more runs would be desirable to 
determine whether different parameter values might lead to very different 
results, Accordingly, four samples were computed, with fifteen random 
vectors each and with combinations of parameters that differed greatly. 
The average number of interview for each RS was given the value of the 
quota of that RS divided by 365 for two runs and twice that much for 
two others, The probability of interview success was fixed by 55% for 
all stations for two runs and 65% for the other two runs. These were, 
of course, combined in such a way as to get ail four combinations. 

Finally, our simulation also produced sample mean loss, an approx= 
imation to the risk for the given estimator, and the sample deviation 
of that loss, Examples of these results are included in Table VII-D. 
When low averages and low probabilities (55%) were used, the percent 
of total quota filled never drcopped below 56% and never reached 57%, 
Individual percent of quotas ranged from 39% to 73%. When low averages 
and high probabilities (75%) were used the percent of total quota 
filled ranged between 75% and 77% while individual quotas were filled 
to levels from 55% to 96%, For the combination of high average number 
of arrivals and low probability of enlistments total quota fulfillment 
never fell below 99% while individual station quotas were filled to 
levels as low as 80% for the very few in each sample not fulfilling 
their quotas. The combination of high probability and high average 
gave 100% fulfillment for all fifteen samples. 

The lack of variation in the percent of total quota filled can 
be viewed as a substantiation of the validity of the statement that 
for large samples jee can be approximated by a normal random variable 


independent of the other m random variables. 
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This lack of variation had one other notable effect. The mean 
B vector, B, proved to be a very good estimator; there were, in fact,no 
samples on which the result was non-optimal, It shculd be noted that 
should the use of the transformation P' result in a sample with a non- 
optimal solution, our program will compute the loss on that sample. 
Therefore, not only was Ba good estimator for our selected inputs, but 
it is also obvious that for these inputs, a static model is adequate. 
Of course, use of different quotas and different parameters might cause 
the probability mass functions governing the inputs to "spread" so much 
that non=optimal solutions might result in enough cases to make a dyna- 


mic model desirable. 


5.4 In the probability model developed herein many simplifying assump- 
tions were made which might have an effect on the conclusions. Some 
of them are discussed below. 

1. It was assumed that quotas assigned to RS's were fixed. It 
has been learned that such quotas may in fact be changed 
during the year in order to capitalize on peculiar seasonal 
or regional situations effecting recruitment. 

2. It Was assumed that the number of interviews held by a recruit- 
ing station could be described by a Poisson probability mass 
function, In using the Poisson distribution, one assumes that 
the probability of more than one event occurring in time At 
is of smaller order than At so that as At approaches zero 
the probability of two or more events occurring in At also 
is negligible. In practice recruiters often give talks to 


large groups of high school or college students, This might 
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be said to be more than one interview; however, it is believed 
that the prospective enlistees must still enter a queue of 
some sort before making a commitment to enlist, so that the 
mass lecture may be regarded as one of the factors determining 
OX, the average number of interviews for the jth RS rather 
than as an actual multiple interview. 

3. Probabilistic independence of the occurrence of interviews 
was assumed, Arguments contrary to this assumption have been 
acknowledged above and some defense given. As a further de- 
fense, we note that the amount of personal effect required of 
the recruiter is a factor tending to promote the statistical 
independence of these random variables. 

4. It has been assumed that the parameters of the Poisson and 
binomial distributions used could be determined for each RS, 
They cannot be completely determined, of course, but it is 
likely that enough data can be obtained from RS records to 
arrive at very good estimates, 

Although some of these objections are forceful, none is believed to be 


strong enough to invalidate the usefulness of the results. 
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Conclusions, 


It can be concluded from this investigation that: 


le 


We have been able to devise a method enabling us to simulate 
(at least approximately) the uncertainties caused by the 
probabilistic nature of numbers of enlistments and thus, can 
consider our model to be more than a mere laboratory affair. 
We have been able to develop a rationale for the derivation 
of an estimator and have been able to obtain an approximation 
of the success we may expect in solving our problem relative 
to an ideal solution, given any particular estimator. 

By means of simulation, we have found that the mean input 
vector, B, is an estimator which yields optimal solutions for 
a very large number of samples and for widely varying para- 
meters. We have seen that this indicates the adequacy of a 
static model. 

We have indicated another promising estimator that might be 
developed that could have considerable appeal to a planner 
under certain circumstances. 

We have been able to devise a method allowing us to present 
to the planner an assessment of the cost of certain decisions, 
namely the choice of assigning various quotas to each of the 
two RTC's, the choice of putting a limit on travel time, and 
the choice of an estimator. 

We have found the simple principle underlying the choice of 


routes to be used by the RS's, that of maximum cost DIFFERENCE. 


24, 








Vey la | 


SWATEONS TO RECRUIT TRAENTNGRC ENTERS 


RECRUITING 


ROM 


COST COEFPFICTEN Ts 


GLAKES SCIEGO 


GLAKES SDIEGO 


GLAKES SDI Rea 


>DIEGO 


GLAKLS 


GLAKES SDTESGO 


~s 


ce 


rn 
les, 
<b 


aa inn He 
sat OFM) 
IT «see 
TOMOO 
(9 ONO 


ae e 


=2=MO 
Cie ~--G 
-— © @ ® 
CAIN A 

IMA 


NW — 
STOO 
“iJ e@ @ @ 
Go) rae 
OMe O 


eee CNY 


rr h oo 
wt OO 
<I 20©¢°e 
oO 0° CXC 

PAM CY 


OR-C 

at ojco 

eo 6¢® @ 
OONGS 
Z2oan™ 
C-— 


deer O'S 
rem i h 
<I eee 
we) 
MOA] 


Ch 
LJ he 
~—> ® « @ 
COW 
CY at CIN 
Tt} 


CAN OO 
Soot oi 
Coe @ 
aA ON 
IOMm st 


CUWQ) 
—OMm— 
e @ e 
a= NIG 
<0O0 
JI—oj 


AMCOW 
Ce 
ee 6 

on aN 
Wht) 


ma 
— 
[a = Ea) 


one CY” - 


<I 


OOO 
m=-— oO) 
CJS © @ @ 
ZONA 
INOW 


aw 


ot ouw) 
LLIN OOS 
wd © © @ 
ON ei 

NA 


OM 
mm LY ON 
~— ee 6 @ 
TONWM 
zOrwo 


Zam 


_ 
Omownw 
PT pit 6 Kea 2) 
— © © @ 
oc — 

(NIOc= 


Koh. 

Cen TDs 
© eee 
COM IG I 
IOsh 
Cj 


al 

Loma 

WONEN 
eo e¢ @¢@ 


Orr: 


LACO 
LUN YES 
oe @ 
COujys 
wed ot 
t—-m 
Le 
WU. Cuh- 
eNCuC. a 
ameee 
at On) 
CUD GAY 


OOOCG 

oo Ji On Bi} 

t @ e 
ma OD 
Coury 


* 
je Ee 


OtMN 
SN COCY 

e oe °¢ 
NO st 
UNCON) 


N 


CIWN 
mm CL SD 
<I r- Co 


WO 
MOON 
Uj e« @ @ 
“Oso 
=O OO 
Cc pe 
= 

Or 
MAOINA © 
LiL ee © 
Sa Ou 

AI 


ON 
ORO 

eo 8 @ 
hc? © 
LULA 


= 
LL sO CON 
OmRNW 
eee 
eye fe 


fe cat 


wo RCM 
—Coe 
eee 

Width 

<I h-O> 

pace | 

ar 

OT onal aaa © @ 

Can OWN 
eee 
aeo-N) 
LWAOiIgy 


CYCy@ 
ONE yaa: 
i ee «6 6 
Dh C70 
Gr (NSO 


= 


> me 


wd af aN 
L1G 
tp ©8¢ @ @ 

wn 


AUN ee 


“ad ol 20) 
UICOhK 
“ e@ 6 @ 
mt) oS 
haa Bi dE era) 
De coal meal 
= 
ada) 
CON — 
(J) @ @ @ 
AD Ope | 
Fos 8 ND 


(a GT) 
ioe a 
Th CD 


mWOO 
> IIL 
p— @e@ 6 @ 
wnt DO 
Lud = CNL) 


INN 
DWAIN 
ae @ 6 
<toOo 
MAIN 


LLM am AS 
wad 
ww © @ @ 
om Or 
MLN CO 


oe cc 


COM IN™ 
MENLO SG 
Cfo e¢ @e 
= C Ora 
gg Fg |) 


CF at or ee | 
met UNM 
—/J @ @ @ 
Co Ores 
0.6 Cop 
gre 
fig 
<Cco 
meh FO Ne 
OD eee 
“Loa Le 


aa ES 


Om 

ee @ 
moony 
Ccr-cw 
— — 
WY 
=> Oe 
a 9 hea 
27 es 6 «© 

0K 


6 Sd ess 


Cu cu ~O 
COG 2 

@ e e 

me NTA 
mt ON OS GS 


hm Orgy 
Se eo) 
C) eee 
CN =f 


CJ[—— 


25 


OID 
WICou) «s 
L.j°©@¢ @ 
a nG 
Ud) i 
ie) 


<= 3} 

Sis 
(M e« © @ 
Oem 
ad ONY 


biG Go 
INCl~ 
(_*‘ © @ @ 
SS a 
ot Okie Tee 

—_— 

xd 
awl Sh -- 
Mm COIS 
j— @© @ @ 
wet PA) Fs VCD 
wd SMA 


me 


raIn 
mC 
<I b- Co 


>~N COW 
CUMIN 
4} © @« © 
OO 

eelta 


CMON 
aot Oth) 
<I ee°® 
wd CS 
Mas ato 


mOID 
MWMOM 
Weeswe 
CL heme ot 
Oars 
> ol aol aes 


SOOM 


LID OO) 


a oe 6 
tous 
Be gs 3 a 


Ore 
NOs 
x= ¢ @ e@ 
NEN 
UIC 
ne | coma eae 
fam 
a7 Oe 

(1G Ye 
"S 0@ @ © 
tus Own 
aot uM) 


(oO =F 
=o 
thi © ee @ 
JQ) Ore 
I= OO 


Shed me 


> 

ae MUI ES) 

SULVLOILE 7 

<I ce @ @ 

zoawn 
NC jr 


as 
Sih OO 
— @ © @ 
JWI OO 
Oger 
(fot —_ 
aca | 
lis =r ONCN 
ZNO OS 
= eee 
aI 
=NNo- 


me 
rN 
mcm 
<= CO 


Mm 


CIOD & 


’ 


ar h-O 
COMIN) ih 
— 
vay 
aNODO 
LOChr- 
(a af eo ¢ @ 
hr 
tON 


NIN) O 
OWN 
eay ee ¢ 9 
Ze Coe 


<tOMmm,. 
me 


CL CO 
CNN 
OX © © 
fh O20 
OM O- 


OIE A 
wh- sO 
(9 «© @ @ 
ewsrot 
DANKO 


ce — 


mw) O8 ot 
he QRWN) 
= 8 & @ 
(oO OB oat | 

ONIN 


it eal © 
HOD0 
ie Oy Hei <@ 
ASN Se 
ad =a 
firm Ale- 


INO 
ZR PCs 
— e e e 
JO a @ Vane 
QAOshn) 


AON) 
— Ce 
[en 
ONY SO 
SU © OT al CO) 
Ite 


2a SO 


CI~ Oo - 


ORIjCN = 
NAMe= 


a 
—" 
“wl NZ 
met CY =~) 
Ihe © 


Tow 
ZO DAI 


Oeee' 


FORM 
COP) fh O 
an=r 
— 
Ler 
Wa Mh) 
<—{I © « »@ 


_ ear 


mon 


AWS 
I 
Ue¢ee 
mm 1 CO SO 
arom 


meas 


~ 


a 
- 


Poe 


Oar te 


Pts rw ter Cate g 


Ye ne ee ee 


lias tn aaa bn 








€8°O 


2) 0A 


GL°2Z 
og°t 
6£°L 
il 
hi °2 
90°E 
geen 


19°L 


tret 
wees 
né*t 


NOLONIHSVM 
SINNOTT LS 
JVLLVAS 
OISTINVAS NVS 
GNOWHIIY 

HOTA Wd 
GQNVIiwOd 
HOYNGSLIId 
VIHd130V1IHd 
VHVWO 

ALID VWOHV 1x0 
AYGA MIN 
SNV317¥0 MIN 


VLIOND ONILINYISY WAANVY TIVY3SA0 JO NOTINEIYISIG 


GO°¢ 
BE °E 
26° | 
Sel 


—ber 


0z°9 
9e°1 
0S°2 
Hh °Z 
99° 
€S°2Z 
geet 
98° 


\ 


JTTVIAHSYN 
SI 7VOdV3INNIW 
334NVMIIW 
NOJVW 
JVTIASINON 
S3139NV SO7 
M90U 3VWLIA 
ALID SVSNVX 
AVTIANOSHOVE 
SI TOdVNVIGNI 
NOLSNOH 


L1OuL30° 


SANIOW S3Q 


II 31798Vi 


ce *t 
ee°t 
€¢ * 
Ln? 
c6°C 
49° 1 
GS °t 
tL? 
ohh 
SO. <¢ 
€& *l 
Meck 
eh * 
9S °2 


YIANAG 
Sv1iva 
SNIWN109 
VIGWnt09 


OGNV1ISZA319. 


ILVNNIONID 
OIVIIHD 
OWjISING 

NOLSO@ 
WVHONIWYIG 
JYOWTL VG 
GNVIHSV 


3NDYANDNG IW 


ANVG TV 


26 


ois 7 oa 


Haile 






ae 
7 -_ 





TAS BE ier! 


OWN OTIRK AN SOM AINITOOMRM—OOKOMN 
MFM — DCOONMNDNO Cer KK NMNMAOMOOD 
oeeereeerteeees#eee#%eesF%% © 6 © @ @ 


NOI AR ANON) AD at StU SOM COR Oren CN SU OR 


ttstastotastst tas ssrnMMMNnWH 


COR OM NID FIOM ORTON AINOKRGTOO 
hh MODOC MHOCMOMOAAAAAAAAAOASO 


ff 


COD COMOROS MH RK OOORIR KH ORE ITOH HM 
= CIN FD OR- OM MDE AR MIND OMN Cr SHMO 
S54 (0 epee te © ee ee eel eel ehlULc lh! Uh! lhe Ce 

MA OR OMDODONSTMNACIMKM OOOO OO to 
DMODONOIMnNIMnnCMnna s-zasastsastas 


RYO FT ODOMDADOMKAVA TNO ODORKNM 40h 
LAUD US NIDA OO 0 OOO OOOO A RE EAAE 


TOON D OF DOW Mm eH KOO ME OMCO 
OM SFOUABRWII AN ODDS CEDMODOODDNOOO 
eeoe#e#se*e?oF#e%eeee ®&e¢ee8e88e@6h6UmM8hmUCUMOMmUCUMAMmUC HDMUC MCL OhUlUhhUhHhUCU OHhmhlUlUhhUCUC Ohm 
= OCOR- IN TMN — OO F- OW SMe OD COOP- OMNI Tt 
HMM BWOCo Mas MO awh-k-BR-KRKK-B-B- OO DOO OOO 


OE-AC OMIM FN OPS ODKOKNM FTWNOM OR Oe 
PIII IIE rm) SS Ss ar a Sara aa on 


OMRAKINOME OMODOD SNK ROK CNM t00R 
OO NKR OMRONRKK AM tOON-KRaeChR t—OnM 


° 8 es eee * e 
Ch NM CUO P-N TN RH OE WY tN OCI O TMA 


QUINN OAIN Ne mm me OO OOS OOOSCAAACGH 


ro.gu on" 9- a2 eww me ee re er Re ee pe ee Pe Tn ee pe 


M— CIM DWN OR-OPAOMAIM TFINOR-DOAOKCIM SW 
CQOO9D ODO 00 D0 cr re rn NNN 


Ott TNOONWIOAONO mK Mr NOMRK OO AA 
—ORDEANOCO FOA ONO CHrMNAMA 
* e@0e@6h6UmemmUCUC“CMCUSGhmUCC HhmhUCUc HhUCc OhUCUC HhUC HMhC TMC HhC<C DMC MWC<‘<‘ HODHCU OULU HhUhUOhUCc OHhlUwF 
FT NOM AON a TING ORE OO OCORmMNM aH 
DONG SOM TIA AON TEOMA MINIM NM) sp op sr or a ot 


4 


CSA OT AINA IL OM ORO MAIN AUTSOR- COO O 
NF-M DOT DO DDOMAARAAAAAHKRAAHO 


t omen 


ee OCC Oe I = OO OO OI OR EO 
MF OR-ARKMODONNM OMVOOANKHMOtOOt 
_e @ @ e808 © @ 8@ © © @ 8 UHmUMOMmhUhCUCUh HhUCUCcOhrhUCc OhmhUCUCUhHhCUCUchOhUCUc HhUCUh hUCUch OhUC OFhUCUCU HhUC Oh 
Mm ODER OW st tlt COOOKP OONN FS Tt 
WO at St TSS tS TFT TON MOO MANTA 


WON O-M EC PRORKAIM AN CR DDORKNM. Ft NOM 
ILD LOG LLL SO 0 000 0 OOO RR EARMR- 


OTOD AANA FOIA Or Or ONE) AOS AAO 
OODDCOCO TK VNANINIM SAS FUIUOO- OO O-— 


OO SNK OOOH ER OIDTMIN H-ORGE™ ON TIMN ° 


RRA E OOD OOO DODOOMNINMNNMNLAW 


OM CRAORKMUIMFNOE DOOCH-NMANOR GAO 
NACIANN HAMANN 4 Ft op ot ott pt tO 


COS SF STUNCIC FNC 3 OOD NOWDOWSMOIN 
OOO NT TOMI MAIN CIN mem mm OO OOOCOCOO 
e* « © @# 8s 6e &@& © 8s @ 86@8ehUrOhmUhHMmUC OMhUCUCc OmhUC<C WC MUlUhTCmUhUhHhUCUCchOl HhUlUhHmhlUhHmhlhUhF 
Ne ORD CHR OWN TMCIMORTA- OM SMO CO OCR: 
ODODARDARAAADASOCWAIMWOCamnmarh-Pr- 


Cr - -_ 


Oe OHH AWD OM OAOMAUM FINCK OC Or CIM Tw 
SOC O22 OOO Ore re ee — NONIOINNE 


27 





oe iY 


Omar em +tWMOIST OR MOMMA NOS Cre 
m= OR OLIN OWCIOR- UM ONES OLY = 
Ce ee eS 2 
SF NPA NIG = SP UV-OOB- OF ODOR AMTUIO 
CUEIN OS CI ON NI ONIN OI OT ION ON QUID NEGRI EOE 


DOOKAINAAOR GCOCRHAIM TNO ORO 
KMS WODODDDDMDMOMAROROAAACRKO 


rc 


COMO St PVO OOM DRONA INNOOWNUH OCDE 
ORE MDOP AIM FOE KDOWNTOOAONIONOOOM 
® e@ee#t#teeteeeee¢ee8e08e8e8 #8 ® @ @® @ @ ® 
MII OP CODON THAOICIR— ODOC OR-CVONS Fst 
Tat SFT MTOR MD MINION INIIAQIONCINN AANA 


NM TOO ODOKNM TNO OO OMNIO TU OM 
LAU LAUI UD LUI 0 0 0 0.0 0 0-0 -QOR ER RAEN 


Woe Oc Osa - OOS ANS Mey OrmnanNk lM 
OCG OUND Ft TNO MOP NAGI st a TUN GOO 


ODOM DATFIGNRHOCOCH OW AMIN K ORGY OW 


MODOC DOODDOOUMNNINUIMON NU St Sot 


OR- BA OM AM SUIOKR- ORO NMANOR- GOO 
oN TaN TaN (0° 10.010 010.0106 1301001).010.01)6 1h 2 Bm aie dike ga oo a an Oe ee 


MN eA CONOR ON FT ST MNINO CN CIN oor QI CINDY 
JIOCCOt as a etc ee 

® ®e ° e @ e e @ e@ 
CONG SNK COON Cin tC OF TOON Cina MN 
OORDRAACHKOOOGOMmUWdocInrmrkR-EEEREES 


OM CIN) IF UNIOP- CID OMAN TUN OM AO Orn ONO FW 
COO OOC DO ee er nt (NCI CIC 


INH OOMMNIOR ANKODK—ORNNOD 
COR-OWAMNOEMOOARCTAMONNRSOMD 
, 
NANA NANA ANON) = Ft TUN CON COROO 
NOUNNNAUNANASN ANANSI 


DMD Om NF a W OR- CID OM NM TUNOR-OAO 
MB-OO CA DOOD CMODAAATAAGKAARANO 


ror 


AOR DOMR ADRK RM OOO OALSR ELINOR Re 
AIR MONIC ARNODRKKRWNOTONRMOAINCID TNO 
e e e e e e e e e e e@ e e e e e e @ t e e L ] @- es e 
MII OR MR OOR CAME MEO OCIA St ot tear 
MIMS FONIMIOAS ENC IOI OI ONION ON OIN AIO AIOIOIOIA 


CUD att -OR- Ss BORD TINO ODOM ANM AUNOM 
UNIDAD NU UL 0006000 050-OR ERR RARR ER 


OO tOODDRE- IM = Or tI NOO TPA FOO TNO 
CAOMUA AQHA CHOSE OLR ARENA Ore AO = OSI COMC 
oe e@® @ @8 @8&@®@6Uh3]PHWDUhC<C<C HWmhC<CwrMhTmChC OOUUCUrhOmhUC HOMUCDMMhCéCrHhC<—CM HFShC<Ci PWOCUrPhmhmCUCh OHUhUCCrOFmhCUCcCOhmC OhUCt HhC H!hC(<C HHSC 
COM IY OW TNO OI OOOO OOM: OONNST tM) 
aes SS Sat ae St SS SIMONI OMIM MIMI) 


OE CIE AORM NIM FAWNOM OS Orr NMA OM OAC 
NANA MSP OM OMANI ANIM Ft St FT TITTY 


UI COON BAIA SH RK AN OM OIMOAOIN tt ~O0 
Passe Sane eae seca COO ON ON 
i) td] t = @ 6 td] e ° . 
cD SNOW SMNHKHRKOCONKE Ou t sre i-OOG 
OOS 06 GOO DOONAN NNN ALUN UMAY) = ot 


Crm OU) STUN CBO OOM NIN FNOR- OAD ORM ANM ALN 
0 OB O90 ODD GO rr ON AIANOINS 


28 








TABLE V 


NUMBER CODE OF RECRUITING STATIONS AND RECRUIT TRAINING CENTERS 
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APPENDIX A 


BASIC SYMBOLOGY OF THE 
TRANSPORTATION PROBLEM 


It is desired to ship a homogeneous product in integral amounts, 
a; (i = 1,...m) from each of m origins and to have it arrive at n 
destinations in such a manner that each destination receives an amount, 
bs (j -1,...n). The cost Cyj of sending a unit amount from the jth 
origin to the jth destination is known for all ij combinations, Here 
the total to be sent ae the total demanded, i.e., 


# 


a 


da, = 


m n 
ist 7 jet 
The problem is to select the amounts Xs to be shipped from the 
ith origin to the jth destinations in such a way that the total cost is 
minimized, There can be no shipments from destinations to origins, 


which is to say Xs j =0O for all ij. Symbolically, we may state the 


problem 


letting: 


m 
Z = 3 ), C15 Xyy (A~1) 


where X4 =O for all i and j 


minimize z subject to the restraints 


m 
yy (A-2) 
Mes = be ] = igen n 
j 9 
ist 
n 
». Ken = & esa m 
7 a gece 
jar“ 
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Xa = 0 for all i and j 


Restating the transportation problem as a standard linear program 
in the matrix algebra form 
minimize 2 = cx (A-1') 
subject to the restraint 


where AX = B (iso) 


¢ =le, J, 


a transposed vector of C3; '8 (all vectors here are assumed to be verti- 


cal arrays, hence, any horizontal array is a transpose). 


i = x, | 


TIGOTe so 
0011 © ° © @) 


A = e we ( A=3 ) 


© 


Os .gecenla 


UC a 


the matrix of coefficients of the system of equations: 


» Xess = Ao 4 Se he om 


, 1) 1 
A 

». Daal c 
; 11 4 


Bie 








we delete the equation for by, since any one of the equations may be 


considered redundant, and 


the adjoined vectors of supplies and demands. 

The investigations of sensitivity can be stated in the following 
ways: 
I. Variations in c! 

Assuming the already-found minimum basic feasible solution is 


Xo, Since it remains a minimum basic feasible solution 


if and only if 


2, = C; for all j 
where = 
xX, = // PR ae 
1 
x; = if P. 


and P"' is the inverse of P, the (m+n) x (m+ n) matrix composed of 
cOlumn vectors of A, 


Then it is true that if 


_ 7 = 
(AC) X5- Ac, = Zs = Cc, 


the solution Xo remains minimum, i.e., in investigating the effect of 
cost on the minimum we first solve (A-1) to see if our old solution is 
still good. If not, we must resolve, 
II. Changes in B 

Changes in B, principally in the partition of it involving the 


amounts to be supplied, are of the most interest in this problem since 
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there is some probability that individual recruit stations will not 
meet their quotas. 
Since the altered B, B' may be described B' =B+ A B= Pt, 


changes in the minimum basic feasible solution satisfy 


Ne ae (eS 


al 


Xi,=P P 


O 


and X4 will be a basic feasible solution, if and only if 


xi > O for all x! members of Xt. 


The changed minimum cost (objective function) 2! is 


i <1 i 
25 = C ee 
cad 
= Zo + C AX, 
: =—T 
Hees At%,=c AX, (A-3) 


Since the values of B are subject to statistical variation, the 
planner responsible for budgeting transportation will want to know how 
such variation will affect his cost, and equation (A-3) gives an indi- 
cation of how to go about it, as we observe later in Chapter 4, 
Statistical Variation. For a complete discussion of sensitivity analysis, 


see Garvin [3]. 
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APPENDIX B 
PROBABILISTIC FORMULATION 

In this discussion, extensive use will be made of the symbols and 
ideas presented in Appendix A and some new ones will be introduced. A 
full list of symbol definitions is presented in the Glossary, Appendix 
Be 

In the problem under consideration we are aided in our probability 
investigations by the fact that of the system of equations describing 
the transportation problem, there is always one redundant equation 
which may be discarded. Since the inputs from the RS's are the members 
of B of most interest, we discard an equation which describes the input 
to one of the RTC's. We now assume that the distributions of the RS 
outputs, Bg)» are independent and of the same form although they 
have different parameters. Of course, the assumption of independence 
may be attacked, as it implies that the events being assigned the 
numerical values b, i.e., the number of enlistments at the RS's, do not 
depend on any common phenomenon such as prosperity, war fever, etc. 
(We do not include the effect of seasonal fluctuations since, although 
they probably affect the rate of enlistment over time intervals shorter 
than a year, they probably have no effect on the annual rate). 

A complication of notation is that when we formulate the problem 


in the simplex notation, we note that 
Be 
Rik 
B= 8; corresponds to} 4a,,| of the 


bj 


transportation notation and the bs of the latter are restricted so that 


das — >, 
4 j 
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or in the notation of our probabilistic simplex 


m m+n 
2. 8:- 2 8 
= j=mt1 


Fortunately in this problem (as indicated above) n = 2 and since the 
(m + ee equation is redundant and is discarded, there is only one 
random variable i, namely S445 that is a function of the others. 


The vector B becomes 





Since we know the proportion ve of the total number of recruits dB, 


which the first Recruit Training Center is assigned, we have written 


ca is 72, B. ; 


In order to obtain the form of the probability distribution func] 
GasOn) sOr B, j = 1,..-.m the following reasoning was used: 

A recruiter has some probability, p, of enlisting any candidate 
he interviews. If these probabilities are assumed identical for any 
one Recruiting Station, the number enlisted in N interviews is a ran- 


dom variable described by a binomial distribution with parameters N 


and p. However, the number of interviews held is largely out of the 
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control of the station, that is, it is controlled in only a general 
way and is not predictable. Energetic RS personnel can raise the local 
level of interest but must wait for candidates to show that interest, 
If time periods are properly chosen it may be said that the number of 
people arriving to be interviewed at an RS is a random variable obey- 
ing the following restrictions: 

1. The probability of one arrival in a period At in length is 
roughly proportional to the average rate of arrivals and to 
the length of At. 

2. The probability of zero arrivals ina period At long is 
roughly one minus the probability for one arrival. 

3. The probability of two or more arrivals ina pericd At is 
of a much smaller order than those in 1 and 2, 

4. The arrival of any one person is independent of that of any 
other. 

These are the standard assumptions for the establishing of a 
Poisson probability distribution law, and the steps used in arriving 
at it are discussed in elementary probability texts (vide Parzen [5]). 
We will simply say that N (the number of interviews held in a RS in 


some convenient time interval At) is described by the mass function 


n 
=t t 
P [N =n] =@ ih sot) n= 0, loose 
= 0 otherwise 


where CX is the average number of arrivals per period (e.g., per 
recruiting day). 
Now the number of men actually enlisted in N interviews is seen 


to be a random variable described by a conditional distribution function. 
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For the jth station 


“Bs (5 [.5) = a = bj = | 


There is one more complicating detail, which is that the re- 


cruiter is not able to stop recruiting until he has filled his quota, 


je In effect, this truncates the binomial distribution governing the 


"output" of his ns interviews, so that 


pe if ny = 93 
and 


Dropping the j subscripts for clarity, then 


FEIN ly) = yy p(b|n) 


{b| p(b| n)>9} 


n 
mi iin? 
K, ye (3) pei 
i=0 
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K, ». (2 | pirn- + 
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+ ey? ae 
n=q+1 
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OS 
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= » D (xt)? ot y Se 
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i=O | n=0 


» (2) pk,.a-k 


k=0 


The size and complexity of this expression suggests the desirability 
of resorting to a computer when computations must be carried out, and 
for some of our purposes we can do so by use of Monte Carlo simulation 
techniques. 

The planner responsible for budgeting the transportation of re- 
cruits will wish to know not only the "best" way to transport them, but 
also what the effort is going to cost. A common measure of cost used 
in such situations as this, where there is uncertainty involved, is to 
determine the expected value of the cost; in this case it would be of 
minimum cost. Ideally the entire probability distribution of minimum 
cost is desired, but problems of this nature are often found to be 
too difficult to solve entirely. In this case even if we could solve 
for the expected value of minimum cost, we would not have arrived at 
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a very satisfactory answer to the question, "How cheaply can this 
transportation be carried out?" The reason the answer is unsatisfac- 
tory is that implicit in the use of the expected value of minimum 

cost to answer the question, is the assumption that a real system 

could be set up to handle any transportation input optimally. In- 
stead of taking such an ideallistic approach we have decided to be 
satisfied with trying to find a solution which will deviate least from 
this ideal, The "analytical” expression for the probability distribu- 
tion of recruiting station inputs to the transportation problem (de= 
rived above) gives some idea of the complexity of an analytical approach 
to finding this expected value. In order to avoid the difficulty, the 
FORTRAN program for solving transportation problems was modified to 
solve a large number of ramdomized problems when given inputs of quotas, 
probability of interview success, and average rate of arrivals for each 
RS, Thus, it was possible to generate a vector of random inputs and to 
obtain a sample mean value of minimum costs. In order to give an in- 
dication of how the distribution of minimum costs varies the program 
was designed to solve for sample deviation also. Here again it should 
be noted that since many data were picked at random (i.e., the values 
of Py the probability of interview success, and (X, the average arrival 
rate for the j*® RS), the solution is only an illustration of method. 


From the previously mentioned relation that 
B|29(B)| =z, (E[B]) (4.1) 


we see that a solution of the transportation problem using mean values 
of inputs will give us an upper bound on the set of possible means of 


optimal solutions, It is also one specific description of the input 
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situation for which we can find a minimum cost transportation pattern, 
but it may not be the "best", 

What we seek for the "best" is first, a loss function describing 
the penalty we pay for improper action, and second, a rule for deciding 
what action to take so that we can expect to pay the least penalty. 

In other words, we are trying to solve a standard statistical decision 
problem and are trying to determine an estimator, Desirable proper- 
ties of such estimators are discussed in [6] and other standard sta- 
tistics texts. Those properties sometimes conflict and are often 
chosen somewhat arbitrarily. Our estimator will be chosen arbitrarily 
also. 


Where we first define 
@ = min 2(B) 
for 8, an estimator of B, a natural loss function,,., would seem to be 


A(B;B) = 2(B;8) - 2, B 


where we write 2(B3B) to emphasize the parametric role played by B. 
later it will also be written without the B. 


The risk, or expected loss, we wish to minimize is 
A 
R(B;B) = £ [A(BsB)| 


Van 
Where B represents the decision rule for arriving at the final action, 


the appropriate transformation on B, 


“n 
Our rule, "use the estimator B", is that which consists of acting 


as though B were the actual value of B, The action then is 
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N> 


N> 
| 
©) 
ae) 
if 
wo 


where we now regard 2 as a transformation on 5, The vector X, the 
"minimum basic feasible solution" is one of the results of the solution 
of our transportation problem under the restrictions represented by B, 
Equation (4.1) becomes 
min R = min E | 2(B) - Zo (B)| 
B 

In words, equation (71) regards the loss as the difference between 
a fixed solution and an ideal solution. The fixed solution is that 
arrived at by use of the estimator 8, The ideal solution is that arrived 
at by use of the actual value of B. 

Drawing on the notation of simplex algorithm one is tempted to say 

min R= min {ore (e[B]) -£ [ot#"5]} 
B 

making use of the linearity of the transformation represented by Ea 
that transformation arrived at by a simplex solution using the estima- 
tor B. In fact, the authors neae this approach until certain computer 
outputs brought the inherent error (discussed below) to their attention. 
Continuing the above reasoning would have led to the following satisfy- 
ing solution, 


letting E[B] = B, and substituting, equation (a) becomes 
g [clp7'p] = c! Bop 
where am is the transformation arrived at by use of B(Blin the simplex 
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algorithm, and by the definition of 


min {2 (B)} = = min che Tale 


aa 


Zz. (B) 


clp-18 s cl p13 


for any other P7! ¢£ Pl, therefore, 


kee eerste Gua: a3 | cp tel} 


A 


B 

ail = {ctr 'B =f [cts] 
Ad = 
Ceo B°= B 
and al = pe 


bf 


the best policy is to use the mean value of B as an estimator and the 
transportation pattern computed from it. 

The above approach is in error in that we assumed we could use the 
transformation 2 | over the entire space of B vectors with some of the 
solutions becoming non-optimal, ignoring the possibility that they 
might become infeasible. However, the fact is that the solution to 
the linear programming problem is either optimal or else it is infeasible, 
where the term feasibility is technically defined. For a solution to 


be feasible it must satisfy the restrictions 


AX 


i 
ee 


IV 
Oo 


and Xs for all i. (Ze) 


For the mathematical details of this aspect of the simplex algorithm, 


the reader is referred to Garvin ily 
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We can restate the definition of feasibility ty defining the 


1 


elements of the P matrix as Pik and writing for equation (4.2) (in 


the case where we started with B and its implicit transformation B71) 


m 
An 
>. Ps a1 
= ik /~k ’ 
and we want to know the probability of this occurring, i.e., 


m 
: 2. vit B.=° 

The implication of the fact that ar EP is optimal or else not 
feasible is that there are two regions in the space of B vectors and 
they are disjoint. In one region the loss function is zero-valued 
and al is optimal. In the other, the loss function takes on a posi- 
tive value, so apparently the estimator B that we seek is that one 
which is associated with the region of infeasibility with the least 


expectation, That is, the risk can be more accurately defined as 


where K is the set of all feasible solutions to the equation AX = B, 
Thus, the estimator we seek is that which minimizes the second term on 


the-rignt, 1.6, 


min R(B;B) = min [{ AGspar, 


P'p gf xy 


A 
B (B|P 
Now, defining X* = ee). and x* as some vector value, 
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Xr. (BB) = 0 if X*E K 


>0 Hepes Kerk 6 mene aK 


<u 
at 
Ps 
K 
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a 
| 
[ 
= 
Ps 
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where K is the compliment of K in the space of all ae). 











and 
fy x (x*) — x, X,, ce has oe x) 
=f (b. .b b )[a(b eon 
ina b 7 gees 9 eee yeood 
jr ek. Lb, ea : 
Now 
J(b,,...b ) = aS diag Ox is Je ¥ 0 
0.» Ob, 
OO an 
db, Ov 
m 
therefore 
m 
a ere at) Ie oe ie Dee eat _ peek: Ie 
Xx Ban B wn = eae? 524 Wa aa: 








Here a difficulty arises due to the non—independence of ome and 


Pe ete ct 
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However, assuming the foregoing operations have been carried out, 
we can make use of the central limit theorem to say that the random 
variables [Gre i = 1,...m, can be said to be approximately normally 
distributed since they are large sums of random variables, at least 
over long periods of time. Further, since 5, is also a large sum 
of these sums, it can be approximated also as a normally distributed 
random variable which is independent of the other m, Consequently, we 


are able to write for these approximations 


fp (b) = ne rene 


iT 


. S ey ™ | 


i.e., there is a density function which can serve our purpose. 


Now 


*IX*E K 


[ A(858)) “dase X) ;B(X*) )dFyx . 


To carry this out we must a % to conform to reality, e.g., 


: Rie - 
Z(B) = & ix? 
' 
where x= 0 if x#<0 
at i 
= x* + x* if x* <0 and i is odd 
a 1=1 
a eee if x* <O and i is even 
at it+1 i+1 
= iy otherwise 


aL 


1 
(This definition of x} models the action of an RS in carrying out a 
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fixed policy). Finally 


m m m m 
Mak (wif ope! iy 
z| Pf sof Ces S/O Seah Da or aie oar fis 
[ i=] alias ee slimes (8 JS ee, = pies iz, mii 
—CO- -CO 


dX 506 dx, 








Here again we have come to a point where for computational purposes, 
theory yields to simulation and the theory merely provides a conceptual 
framework for what we are trying todo. Of the ideas written above, at 
least the loss function can be useful in the simulation, 

In investigating the results of the above reasoning by Monte Carlo 
Simulation, a program was designed to provide the randomized input B 
vectors and solve for their sample mean (one common estimator). This 
mean vector was then used as an input. The program computed the results 
of transforming the randomized B vectors (using the transformation found 
in the solution for the mean value vector) altering it by using the de- 
finition given above for xe outside the region of basic feasibility. 

It further computed the loss function in each case and the sample mean 
loss or sample risk for this estimator, In computing the random B 
vectors, the same concept of the events giving rise to enlistments was 
used as in the analytical description; however, because of the large 
numbers involved, and to speed computation time, a normal approximation 
to the Poisson distribution was used. For this approximation see 

Parzen [5], and for the formula used to compute random numbers according 
to a normal distribution see Vaa }9|. 

We are able to find sample risks given specific B estimators but 
we do not believe it is possible to solve for the best. The only one 


for which we have solved so far is the mean B vector, B. For examples 


oS, 








of the results of this simulation (using completely arbitrary para- 
meters) see Tables VII-A through VII-G, pages 55 to 61, and for a 


discussion of the results see Section 6. 
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APPENDIX C 

OTHER POSSIBLE LINES OF INVESTIGATION 
1. When a planner decides to seek a solution to a problem by using 
the criterion of minimum expected loss, he is considering all events 
(in our case B vector values) which have any positive probability of 
occurring. Suppose, however, he should prefer to make a wager based 
on the assumption that in the finite period of time his system (the 
U. S. Navy) is in operation, that events of extremely small probability 
will not occur, By basing his actions on such an assumption he is bet- 
ting against nature that he can save an amount of money in this finite 
period of time, realizing that he risks losing more than the amount he 
saves, but risks it with a very low probability that the (extra) loss 
will occur. In effect, in our case he wishes to find the solution 


(B estimator) for the following objective function: 


(38) A } 
i | 2 (B3B) = 2,(B)| dF 
"~ Lee Zo L 


where 
Q represents the upper limit of the small level of probability of 
the events he wishes to bet will not occur, 
% (B;B) is the objective function solution of the transportation 
problem for B evaluated at B and emphasizes the parametric role 
played by 8. 
dF is the differential of the probability distribution function 
for A(B;8). N(B;B) 9 is the Ot percentile of A(B;8) under the 
probability law F,. 


That the bet mentioned above might be a desirable one to make 
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could be determined in the following way. Suppose, as is likely, that 
the expected values of the B's, the components of our random vector B, 
occur near their maximum values, the q's. The reason this is likely 
is that the quotas of recruiting stations would probably be assigned in 
just such a manner, Also suppose B is our best estimator using expec- 
tation. Further, suppose a sensitivity analysis showed that small 
changes of B With values of £ increasing should cause the solution to 
become non-optimal, i.e., the loss for small increments of (5 increases 
rapidly and should not satisfy the constraint that x; = 0 for ail x, 
of the solution; and that equal small changes of B with values of 5 
decreasing should not cause the soiution to become non~optimal, i.e., 
for small decrements of §, the loss remains at zero. The probability 
mass function of one of such B's, ihe would have a form similar to 


Figure A, 


(b) 


fe 


FIGURE A 


From the figure we see that the probability that b takes on 


the exact value (5 + OB) is greater than the probability that it 


takes on the exact value G -A fs) 
BB * 48)" % 8 -4p? 
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where A Bis some arbitrary small change in B, This statement hinges 
on our assumption of the level ot B relative toq and the fact that 
q is large, i.e., that B is on the increasing side of the mode of 
p,(b) and so are both (B +O) and (B -Af). 

Looking in the same general way at the loss function, we see that 


it will probably have a form similar to that of Figure B, 


A(B3B) 


Bb * B 
FIGURE B 


since there will obviously be no loss an Ble Tt might (and a sensitivity 


analysis would determine this) Look like Figure C, 


A(B3B) 


fe tp 
FIGURE C 
so that, in the latter case the probability density function of \(B;8) 


might be similar to Figure D. 
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( \ (B;B)) 


Bp «a B 
FIGURE D 
Even so simple a discussion as that above makes it clear that, 
should these conditions hold, for all (or possibly even some smaller 
number) of the components of B, an estimator of B with components to 
the right of the prospective Bis would be a reasonable choice, given 
the willingness to make the previously mentioned bet against nature, 


The choice of such a new estimator, B, might result in a loss function 


such as that in Figure E, 


A(B3B) 


ms p 


and a probability mass function similar to that of Figure F, 
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p( A(BsB)) 


BR 2 fp 
FIGURE F 

The small standard deviations resulting from the parameters we 
used in our simulations, indicate that this second approach would not 
be profitable, However, were the numbers of recruits sufficiently 
small to "spread out" the probability mass functions more, an inves- 
tigation of this additional appreach might become profitabie. 

The effect of ignoring the tail, that is, those losses with 
probability less than 6, is illustrated in Figure G. In effect the 
new B, B, can be found by renormalizing BR, truncated at the eee 
tile and finding the new mean of the region of B vector which produce 
the truncated function. An approximation to this new mean (a new 


sample mean) could easily be found by use of a computer, 





FIGURE G 
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2. A possible improvement that might be made in the solution of 

this problem of the transportation probability model would be obtained 
by investigating the problem's dynamic aspect. A more sophisticated 
model than ours might allow the planner to make an optimal decision at 
any point in the time period (one, n, or an “infinite number of years) 
based upon the conditions then prevailing, It is believed that a multi- 
staged linear program such as is discussed by Dantzig [4 | and [7 | or 

a dynamic program approach such as those discussed by Bellman [8 | could 


be used to provide this dynamic aspect, 
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APPENDIX D 


GLOSSARY OF SYMBOLS 


The matrix of coefficients of the system of linear equations 
represented by equations (A=2) and (A=2'), The matrix is 
defined in (A-3). 


The total number of recruits to be shipped from the jth RS, 
of which there are m, 


The average number of candidates arriving per day at the 
je tot 


The adjoint column vector composed of the a, 'S and Dee 

A fixed value estimator of B, 

The vector of expected values of the random variables (Bie 
A (second) fixed value estimator of B. 


The components of the vector B when the vector is considered 
as a vector of random variabies, 


A fixed value in the domain of bie Aliso the total number 
of recruits to be shipped to the yeh RTC, of which there 
are gn. 


The transposed column vector of See 


The total cost of shipping one recruit from the 4th RS to 
the j“ RTC, 


The proportion of total input assigned to the RTC of interest, 


The expected value of some random variabie V in which V might 
be a scalar or vector, In the iatter case, the symbol re= 
presents a vector of expected values. 


The probability distribution funetion of some random 
variable V. 


The Jacobian determinant of a transformation with respect to 
the variables (...). 


The “a of all feasible solutions to the linear program for 
1 Sa 3 


1 


The complement of K in the space of all possible jor as 


vectors. 
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a 


ms 
eke Che 


ra 


A loss function, the difference between the least cost 
possible if the problem could be solved optimally for B, 
and the cost obtained using the transformation derived by 
using B as the estimator, 


The number of interviews occurring at the jth RS. 
The inverse of the matrix composed of a set of linearly 
independent vectors, used as a tran:formation on B to 


yield X. 


Those transformations arrived at by use cf BR and B respece 
tively in simplex solutions of the transportation problem, 


The ae element of matrix P', 
The probability mass function of some random variable V. 


The probability that the jth recruiter will enlist a candi- 
date at an interview. 


The vector of qj'S- 

The assigned quota of recruits for the us no: 
A regret function or expected value of loss. 
ee i Oe 


J J 


The level of probability of the events which the planner 
bets will not occur. 


The column vector composed of the Xq 7/8 


The total number of recruits shipped from the ith RS to 
the j&® rtc, 


The vector of x;4's which yields the minimum z for the 
general case. 
An element of X,. 


a 
The vector of x,4's which yields the minimum z when B is 
the vector of restrictions, 


The vector valued random,variable resulting from the operation 
si | 
of the transformation on the random vector B, 


A possible value of X*, 
th A 
The i” component of X*, 


The 1%h component of x*, 
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N> 


N> 


iit 4 


The transformed x# which models the actual situation of an 
RS carrying out a policy. 


The objective function expressed in equations (A-1) and 
(A-1'), 


The transformation which yields the minimum objective 
PuUNCE TC ees o eee. 


The minimizing transformation found using B as a get of 
parameters and considering B as a vector of random variables, 


The minimizing transformation found by using fo as the speci=- 
fie value of B, 


The transformation of B which gives the minimum of the 
objective function when acting on B. 
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